Testing Lorentz Invariance Violation with the WMAP Five Year Data 
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We consider different renormalizable models of Lorentz invariance violation. We show that the 
limits on birefringence of the propagation of cosmic microwave background photons from the five 
year data of the Wilkinson microwave anisotropy probe (WMAP) can be translated into a limit of 
Lorentz symmetry violation. The obtained limits on Lorentz invariance violation are stronger than 
other published limits. We also cast them in terms of limits on a birefringent effective photo "mass" 
and on a polarization dependence of the speed of light. 
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I. INTRODUCTION 

The principal spacetime symmetry of particle interac- 
tions in the standard model is Lorentz invariance. Ex- 
periments confirm Lorentz symmetry at all currently ac- 
cessible energy scales of up to 2 TeV. This scale will 
be extended shortly to 14 TeV with the Large Hadron 
Collider (LHC) at CERN. Although present experiments 
confirm Lorenz invariance to a good precision, it can 
be broken in the very early Universe when energies ap- 
proach the Planck scale. There are a number of exten- 
sions of the standard model of particle physics and cos- 
mology that violate Lorentz invariance (for reviews see 
Refs. |1 1 1 Hi). 

As it can be expected, Lorentz invariance violation 
(LV) affects photon propagation (the dispersion rela- 
tion), and generically results in a rotation of linear po- 
larization (birefringence) . Other effects include new par- 
ticle interactions such as a photon decay and vacuum 
Cherenkov radiation [J|. All these effects can be used 
to probe Lorentz invariance. The dispersion measure 
test is based on a phenomenological energy dependence 
of the photon velocity [6| (see also Refs. pj for reviews 
and Refs. @, S [lO| for recent studies of this effect; early 
discussions include Refs. [1JJ; Refs. p|, @, LL1| consider 
Lorentz symmetry violating models which preserve rota- 
tional and translational invariance but break boost in- 
variance). 

Several models of LV predict frequency dependent ef- 
fects. Such high e nergy L orentz invariance breaking are 
discussed in Refs. [12|, [Dl, LLJ] ■ Refs. [15|] study general- 
izations of electromagnetism, motivated by this kind of 
Lorentz invariance violation. On the other hand, LV as- 
sociated with a Chern-Simons interaction [16|, [TtJ affects 
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the entire spectrum of electromagnetic radiation, not just 
the high frequency part, and induces a frequency inde- 
pendent rotation of polarization (see Sec. 4 of Ref. [2(] 
and Sec. Ill of this work). 

To determine the effects induced by Lorentz symme- 
try violation, it is useful to consider the analogy with the 
propagation of electromagnetic waves in a magnetized 
plasma as outlined in Refs. @, QJl, [Ty, [l8|, [jj|. Using 
the well known formalism for the propagation of light 
in a magnetized plasma, is easy to see that for Lorentz 
symmetry violating models which depend also on polar- 
ization and not only on frequency, the rotation measure 
constrains the symmetry breaking scale more tightly than 
the dispersion measure, see Refs. [14, Ha. Il9l] . 

The propagation of ultra-high energy photons repre- 
sents a promising possibility to probe Lorentz symme- 
try [20| . Gamma Ray Bursts (GRB) are astrophysi- 
cal objects located at cosmological distances which emit 
very energetic photons |6[; reviews describing cosmologi- 
cal tests involving GRBs are e.g. Refs. 0,[21(, for recent 
studies see [22j . After the observation of highly linearly- 
polarized 7-rays from GRB021206 has been reported [25| . 
Refs. [2J, [2J| have proposed to test Lorentz symmetry 
violation with the rotation measure by analysis of GRB 
polarization. Even though this measurement has been 
strongly contested [26J, there is evidence that the 7-ray 
flux from GRB930131 and GRB960924 is consistent with 
more than 35% and 50% polarization, respectively [27j |. 
However, the issue of polarization of GRB 7-rays is still 
under debate and additional X-ray studies are needed to 
either confirm or disprove polarization of 7-rays from 
GRB's [1|. 

In this paper we mainly consider renormalizable mod- 
els of LV as described in Ref. [2| ■ We use the very well un- 
derstood and measured temperature anisotropy and po- 
larization of the cosmic microwave background (CMB) to 
constrain Lorentz symmetry violation. These data have 
been proposed as a probe of Lorentz invariance in the 
Universe in Refs. [29|, [30, bill (32j. In our study we use 
the WMAP 5 year limits on birefringence [33| and obtain 



limits which are significantly more stringent than those 
obtained from radio galaxy polarimetry [16]. 

As we shall see below, generically Lorentz symmetry 
violation leads to birefringence, i.e. a photon disper- 
sion relation which depends on polarization. This leads 
to a rotation of the CMB polarization which induces 
parity-odd cross correlations, such as Temperature- B- 
polarization and E'-iJ-polarization [29j. These correla- 
tors vanish in models which preserve parity. Generally 
speaking, the effect is similar to that induced by a ho- 
mogeneous magnetic field [33, |35[ . In this paper we use 
the WMAP-5 year limit on the rotation measure [33| to 
contrain Lorentz invariance violating theories. 



II. LORENTZ INVARIANCE VIOLATION: 
GENERAL DESCRIPTION 

For methodological purpose let us first briefly summa- 
rize the usual Faraday rotation effect. We consider an 
electromagnetic wave with frequency u and spatial wave 
vector k, k = |k| propagating in a magnetized plasma. 
A linearly polarized wave can be expressed as superposi- 
tion of left (— ) and right (+) circularly polarized waves. 
In a magnetized plasma, a homogeneous magnetic field 
induces a difference in the phase velocity of left and right 
handed waves. This causes a rotation of the polariza- 
tion, called Faraday rotation [3a ]. The group velocity of 
the wave also differs from c. These two effects can be 
expressed in terms of the refractive indices defined by 
k± = n±uj where k± denotes the wave number for right 
and left handed waves. The indices n± are 361 ] 
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Here u> p = 4ire 2 n e /m e is the plasma frequency and 
lu c = eB/m e is the electron cyclotron frequency for the 
magnetic field B (see Sec. 4.9 of Ref. [361]). 

The magnitude of both the dispersion measure, due to 
the different group velocities, and the rotation measure, 
i.e., the rotation of polarization, are proportional to the 
photon travel distance A/, 
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Here, At± is the difference between the travel time of a 
right-handed (left-handed) photon and that of a photon 
traveling at the speed of light, and Aa is the rotation of 
the angle of polarization. 

Faraday rotation is widely used in astrophysics to mea- 
sure magnetic fields in galaxies and clusters (see Ref. (37| 
for a review and references therein) . In cosmology, Fara- 
day rotation of CMB photons [3J, [38| has been used to 
constrain the amplitude of a homogeneous as well as a 
stochastic cosmological magnetic field [39j . 



In the following, we show that Lorentz symmetry 
violation leads to a modification of Maxwell's equa- 
tions [lj, LL5I] analogous to the modifications described 
above. 

Following Ref. [4], the most general renormalizable 
form of Lorentz symmetry violation can be expressed by 
two additional terms in the action (we set H = c = 1) 
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where Greek indices (/i, v, A, p) denote time-space coor- 
dinates, F^ is the electromagnetic field strength ten- 



sor, F„ v 



£y,v Xp F\ p is its dual, 
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is the totally an- 



tisymmetric tensor normalized such that £0123 = \J — 9 
and A v is the vector potential. The four-vector (L M ) = 
(Loi L) has the dimension of mass and describes a super- 
renormalizable (dimension 3) coupling and Kp U \ p is a 
renormalizable, dimensionless coupling giving raise to a 
dimension 4 operator. We want to break Lorentz symme- 
try, but keep conformal invariance of electrodynamics in 
this work. For this we have to ask that the components 
of Kp V p and Lp be independent of conformal transfor- 
mations of the metric. In the cosmological setup with 
9nv = a 2 (t)r] l j lU , the above action is then independent of 
the scale factor a(t). I.e. in a conformally flat spacetime 
the action is like in flat space. To see this note that the 
forms Ap and Fp V are independent of the metric hence 
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, jF p.v F x P and L^A v Fp V scale like a 
celed by \/—g — a 4 . 

The tensor Kp V \ p has the same symmetries as the 
Riemann tensor and we only consider its trace-free part 
which is analog to the Weyl tensor (the trace part also 
leads to dispersion measure but not to birefringence, we 
therefore do not consider it here). Even though we apply 
the formalism used for the Weyl tensor below, we do not 
consider Kp V \ p to be the Weyl tensor which of course 
vanishes in a (unperturbed) Friedmann universe. The 
most plausible origin for the Lorentz violating terms in 
([4]) is that Kp V \ p or the vector Lp stem from the non- 
vanishing vacuum expectation value of some dynamical 
field and the action ([3]) therefore represents a sponta- 
neous rather than explicit breaking of Lorentz symme- 
try. However, for the following discussion the origin of 
the Lorentz violating terms is not relevant. 

Both terms in Eq. ^ lead to birefringence but the 
frequency dependence is different. The first term in the 
action Tlv can be computed within the Newman-Penrose 
formalism, which is usually applied for the Weyl ten- 
sor [2]. We consider a plane wave with conformal wave 
vector (kp) — (w,k). We normalize the scale factor to 
unity today, ao = a(to) = 1, so that conformal frequen- 
cies or length scales correspond to physical scales today. 
In terms of the conformal wave- vector, the dispersion re- 
lation is like in flat space where it has been derived in 
Ref. i, 
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Here \l/o is the analog of the Newman-Penrose scalar (for 
more details see [2fl), 



which yields n± = 1 ± r y(u>)/2. Using Eq. Q, we obtain 
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where m and m represent the left and right circular po- 
larization basis vectors and n = k/fc is the photon prop- 
agation direction. We normalize n and m with the flat 
metric, n l n J Sij = rrffh^ij = 1, and multiply the expres- 
sion with the correct power of the scale factor, a -4 , so 
that, given the scaling of the tensor K, one sees explicitly 
that ^o is independent of the scale factor. (Latin indices 
indicate spatial components of a vector or tensor.) 

The second term in the action IYv leads to the disper- 
sion relation 0, LL6J 



{k^f + {k^){L v L v ) - {L^f = 0, 


(6) 


To first order in the small parameters L 


and 


y^S^LiLj = L one has 




J 2 = k 2 qp w(L — L cos </>), 


(7) 



where <$> is an angle between the photon propagation di- 
rection and the vector L, cos</f> = (L • n)/L. Note the 
similarity of the expressions ([5]) and (O with the corre- 
sponding ones following from Eq. (fTJ). 

To be as general as possible, we rewrite the dispersion 
relation for both types of Lorentz symmetry violation in 
the form (see also fl9j|L 
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where M PL is the Planck mass, M PL ~ 1.2 x 10 19 GeV, 
N is the dimension of the Lorentz symmetry violating 
operator and M is a mass scale of the model. For TV = 4, 
the birefringent part is independent of the photon energy 
and we have 8* = M/M PL . For N = 3 the Planck 
mass cancels out and the birefringent term is inversely 
proportional to the photon energy. The mass scale is 
M = Lq — L cos <j>. Generally speaking, the smaller M, 
the weaker LV. For N = 4, LV is frequency independent 
and the amplitude of the effect is of order M/Mpj, while 
for the super-renormalizable case, N = 3, LV is strongest 
at low frequencies, u) < M. Our aim is to limit the 
function 



Y(w) 
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from CMB birefringence. This ansatz can also be applied 
to non-renormalizable models with higher dimension op- 
erators. For N > 5, M ^ indicates that there is LV at 
frequencies w >Mpi(Mpi/M)'* r = 3 . 



III. RESULTS 

To compute the CMB polarization rotation angle in- 
duced by Lorentz symmetry violation, we follow the anal- 
ogy with photon propagation in a magnetized medium 



Act' 11 ') = -w7(w)AI. 



(9) 



In the case N — 4, 7 is frequency independent, hence 
Aa^ LV ' grows linearly with frequency. In this case, and 
for all models with higher dimension operators, the best 
limits can in principle be obtained from high frequency 
photons (for example GRB 7-rays [H, H3]), while CMB 
photons are less affected. However, the fact that the 
theory of CMB anisotropies and polarization yields that 
both TB and EB polarization have to vanish in stan- 
dard cosmology, while the polarization of GRB's is still 
under debate, at present, a test using CMB data is to be 
preferred. Another advantage is that for the CMB the 
distance Al ~ H^ 1 is maximal. 

In the dimension 3 model, Aa^ LV ^ = — ^(Lo — 
Lcos4>)Al, is frequency- independent. In Ref. [16| the 
above result is applied to polarization data from distant 
radio galaxies, Aa < 6° at 95% C.L. at redshift z ~ 0.4. 
The constraint obtained if Ref. [16fl is \Lo — L cos cj)\ < 
1.7 x 10~ 42 ho GeV, where ho — 0.7 is the present Hubble 
parameter in units of 100 km s _1 Mpc . 

We use the recent WMAP-5 year constraints on the 
rotation angle of the CMB polarization plane (combined 
constraints from the low and high multipolc number, 
I, data assuming a constant Aa across the entire mul- 
tipole range), -5.9° < Aa < 2.4° at 95% C.L. and 
Aa = -1.7° ± 2.1° at 68% C.L. [H (Sec. 4.3). Assum- 
ing Gaussian errors, it is straightforward to convert this 
to the following limits on the absolute value of rotation 
angle, 

|Aa| obs < 4.90° at 95% C.L. , (10) 

|Aa| obs < 2.52° at 68% C.L. . (11) 

We adopt Al ~ 9.8 x lO 9 /^ 1 years. We express our 
results in terms of i^ioo = f/100GHz to keep them as 
independent of the CMB band frequency as possible. 

Using Eq. ([5]), we find the following limit on the func- 
tion j(v) with lu = 2ttu: 

7(1/) < 8.6 x 10- 31 ^- o/io at 95% C.L., (12) 
7(1/) < 4.4 x 10 _31 ^/i at 68% C.L., . (13) 

We can also express the limit on 7 in terms of a limit for 
the mass scale M or the dimensionless parameter M/Mp\ 
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Mm 

(14) 
For N > 4, these limits are not very interesting, while 
for N — 4 or N = 3 'naturally expected' values of the 
parameters are ruled out. More precisely, for the models 
considered we constrain the dimensionless scalar ^0 f° r 
the N = 4 case, 

|*o| < 1-1 x 10 _31 fto^r o at 95% C - L -' 



while we find for N = 3 

\L -Lcos4>\ < 3.6 x lO- 43 h Q GeV at 95% C.L. . 

This is almost an order of magnitude better than the 
limit obtained in Ref. [lta ]. 

We can also introduce an effective photon "mass" by 
writing the modified dispersion relation in the form lo 2 = 
k 2 ± m 2 with 



m = uj j(u>) = Mui 
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For N > 2 this is not a mass in the usual sense 
of the energy of the particle at rest, but rather a 
measure for the modification of the dispersion relation 
which tends to zero with frequency. For the renormal- 
izable dimension 4 and 3 operators considered in this 
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work we have 

[u)(\Lo ~ L cos 4>\)Y'^ respectively. As in Ref. [23JJ we 
can interpret our result also in terms of a polarization 
dependent group velocity, 



v± = l± 



N-2 M ( u> 



2 Afpi V Mpi 
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Ref. [23j only studied the cases N > 5. From Eq. (|12p we 
derive the constraint on the effective birefringent mass, 

777 7 < 3.8 x 1(T 19 (h v W0 ) 1/2 eV at 95% C.L. (16) 

Note that left and right handed photons have effective 
square masses of opposite sign. For the velocity difference 
this implies 



\v+-V-\ < 



8.6 x 10- 31 h v 100 
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00 
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100 



at 95% C.L., for N = 3 
1.7 x 10- 3U /i uroo at 95% c - L -> for N = 4 

The limits on 77i 7 are model independent because m 7 only 
depends on the directly measured rotation angle Aa and 
on the frequency. 

If L <C Lo, we can safely neglect the angular depen- 
dence, and assume that m 7 = \[ZjLq. However, if 
L > io, the modification of the photon dispersion be- 
comes direction dependent, and must be averaged over 
all sky for the CMB photons. Then, the rotation angle 
can be estimated by the two-point correlation function, 
i.e., Aa e ff = (lAal 2 ) 1 / 2 . A rough estimate leads to a 
pre-factor ~ l/\/2. In a more detailed analysis the pres- 
ence of L breaks rotational symmetry and leads to off- 
diagonal correlations in the temperature anisotropy and 
polarization spectra analog to the effects on the CMB by 
a constant magnetic field [35, 38]. To take this fully into 
account requires to estimate the CMB Temperature-B 
polarization, E- and B-polarization cross correlations, as 



well as B-polarization spectra due to the Lorentz symme- 
try violating vector field L, and to compare theoretical 
estimates with the corresponding CMB anisotropy and 
polarization data. Also the scalar |\to| of the 4D model 
breaks rotational symmetry and taking the direction de- 
pendence of Aa into account is relatively complicated. 
This breakdown of statistical isotropy can also be tested 
using the bipolar power spectrum introduced in Ref. 45] . 
We shall address this issue in future work, but even 
though the limits may improve somewhat, we do not ex- 
pect them to change significantly. 



IV. CONCLUSIONS 

The obtained bound on a birefringent effective photon 
mass is below the limit for a standard photon mass given 
by the particle data group [4lj], to 7 < 3 x 10 -19 < 10~ 18 
eV, but is less stringent than the limit from galactic mag- 
netic fields which is, however model dependent [42]. Of 
course our photon mass would be measured only when 
measuring the dispersion relation of a polarized photon 
beam and would disappear when averaging over polar- 
izations. It is not an ordinary mass. 

Another useful bound is the departure of the refrac- 
tion index in vacuum from unity, i.e., \ An\ = |1 — k/u>\ = 
|7M|/2. In the 4D model, |Ati( 4D )| ~ 4|* |, In the 3D 
model, |An (3D) | ~ L /2uj (when L < L ). Generically 
Eq. HID implies |An| < 4.3 x lCr 31 ^^. The differ- 
ence of the refractive index from 1 can be viewed as a 
difference of the photon speed from 1, Ac at the level of 
1CP 30 , which is much more st ring ent than the (more gen- 
eral) limit obtained in Ref. [43], which is Ac < 10 -23 . 
The formalism given here is applicable also for higher 
dimension operators, but due to the frequency depen- 
dence |a (Ly) | ex /" 3 , the CMB based limits on the 
amplitudes for higher dimension operators become much 
weaker that those given from high energy photons (7 
or A-rays). Even the bounds obtained from the nearby 
Crab Nebulae are more promising [4J] if N > 5. 
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